Generation of time-bin entangled photon pairs using a quantum-dot cavity system 
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We present a scheme to realize a deterministic solid state source of time-bin entangled photon 
pairs using cavity-assisted stimulated Raman adiabatic passage (STIRAP) in a single quantum dot. 
The quantum dot is embedded inside a semiconductor cavity, and the interaction of a coherent 
superposition of two temporally separated input pulses and the cavity mode leads to a two-photon 
Raman transition, which produces a time-bin entangled photon pair through the biexciton-exciton 
cascade. We show that the entanglement of the generated state can be measured using triple coin- 
cidence detection, and the degree of entanglement is quantified as the visibility of the interference. 
We also discuss the effect of pure dephasing on entanglement of the generated photon pair. Pro- 
nounced interference visibility values of greater than l/y/2 are demonstrated in triple coincidence 
measurement using experimentally achievable parameters, thus demonstrating that the generated 
photons are suitable for applications with Bell's inequality violation and quantum cryptography. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Dv 



I. INTRODUCTION 

A source of entangled photon pairs is an essential build- 
ing block for various quantum information processing 
protocols^, such as quantum cryptography^ and quan- 
tum teleportatioiA Generally, the employed entangled 
state of photons in these experiments are entangled in 
both the energy and the polarization degrees of freedom*. 
However, because of unavoidable polarization dispersion 
in optical fibers, the polarization entangled photons are 
not suitable for distribution over large distances. In re- 
lated experiments^—, entangled states of photons in en- 
ergy and time degrees of freedom, using discrete time in- 
terval (time-bin) for photon emission, have been demon- 
strated and the entanglement between these photons has 
been successfully distributed over distance of 50 km 8 -. In 
these experiments, the time-bin entangled photons were 
generated through the parametric down convertor (PDC) 
using the pump as a superposition of two time separated 
pulses. A PDC is a heralded source of entangled pho- 
tons where the number of generated photon pairs are 
probabilistic^. At low pump intensity, when the proba- 
bility of generating more than one photon pair remains 
small, the efficiency of the source remains very low (less 
than 20%£). For quantum information processing appli- 
cations, one requires a scalable source which generates 
precisely a single photon pair on demand^. In the last 
few years, there has been considerable progress for de- 
veloping on demand single photon and entangled photon 
sources using single quantum dots (QDs)^i~— , where the 
QDs provide the potential advantages of integrability and 
scalability in such experiments. In semiconductor QDs, 
polarization entangled photons have been successfully 
generated in the biexciton-exciton cascade decayi^— . 

In 2005, Simon and Poizal?^ proposed an on-demand 
generation of time-bin entangled photons through the 
biexciton-exciton cascade in idealized QDs, where the 
bi-exciton state is created by pumping through two 
pulses interacting at two distinct times. The state 



of the time-bin entangled photon pair is given by 
IV') = \/Pi\ early) 1 1 early) 2 +e l ° v / p 2 ( 1 - Pi ) | late) 1 1 late) 2 , 
where early and late are two time bins and pi is the prob- 
ability of generating a photon pair in the early time bin 
(from the first pulse) and p 2 is probability of generat- 
ing a photon pair in the late time bin (from the second 
pulse); the total probability is then pi + p 2 = 1- For 
generating maximally entangled state \ip), one requires 
Pi = p 2 = 1/2. Therefore, a precisely regulated popula- 
tion transfer between the QD energy levels is essential. 
Moreover, pure dephasing processes present in semicon- 
ductor produce detrimental effects on the entanglement 
of the generated state. In quantum information proto- 
cols, such as entanglement swapping, it is essential that 
the photons should not have any other correlation except 
the time-bin entanglement. However, in the biexciton- 
exciton cascade, emitted photons also have time corre- 
lations. These undesirable temporal correlations can be 
minimized by manipulating emission rates of photons us- 
ing resonant cavities^. 

In this work, we propose to generate an efficient time- 
bin entangled photon pair using stimulated Raman adia- 
batic passage (STIRAP). The coherent excitation in the 
system of QDs embedded in a semiconductor cavity have 
been an active area of researc h 17 ' 18 . We consider the in- 
titial QD state is in a mctastable state, and there have 
been several methods for achieving this using electrical 
control of QD-cavity mode resonc o 19 ' 20 . We demonstrate 
that the STIRAP process then provides an efficient reg- 
ulated way for population transfer. We also investigate 
how the cavity enhanced decay rates suppress the detri- 
mental effects of pure dephasing. 

Our paper is organized as follows. In Sec. II, we present 
a formal theory of generation of time-bin entangled pho- 
ton pair from a single QD coupled to a semiconductor 
cavity. In Se. Ill, we investigate the measure of photon 
entanglement by a triple coincidence detection and also 
study the effects of dephasing. In section IV, we present 
our conclusions. 
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II. GENERATING TIME-BIN ENTANGLED 
PHOTON PAIRS USING STIRAP 

We consider a QD embedded in a semiconductor mi- 
crocavity, where the energy level diagram of the system is 
shown in Fig. [TJ The dipole transitions from the biexci- 
ton state \u) to the exciton state \y), and from the exciton 
state \y) to the ground state \g), are coupled through a y- 
polarizcd single mode of the semiconductor cavity, with 
coupling constants g% and <? 2 j respectively. Because of the 
large biexciton binding energy of semiconductor QDs, it 
is not possible to couple the biexciton and exciton transi- 
tions from the same cavity mode, and thus manipulation 
of the biexciton binding energy becomes essential in these 
system a 21 ! 22 . Usually the binding energy of the charge- 
neutral biexciton has a negative value, however by chang- 
ing the confinement size 2 ^ or by changing the strain 2 ^, it 
has been found that the biexciton binding energy can be 
tuned to zero or a positive value. Very recently, manipu- 
lation of the binding energy of the biexciton has also been 
reported by applying lateral electric fields^ 2 -. Moreover, 
construction of an electrode for applying a lateral elec- 
tric field in the vicinity of a QD within a photonic crys- 
tal cavity has also been reported^. Therefore, it is now 
possible to manipulate the binding energy of biexcitons 
inside semiconductor cavities. The another advantage of 
lateral electric fields is that they can be used to create a 
voltage-tunable metastable exciton state \m). 

In III-V QDs, bright netral excitons are formed when 
an electron from the s— shell or the p— shell of the valance 
band is excited to the s— shell or p— shell of the conduc- 
tion band, respectively; the transitions from s— shell to 
p— shell and vice versa are essentially symmetrically for- 
bidden. However, in the presence of an applied lateral 
voltage, the charge carrier symmetry can be suitably 
broken broken so that a bright exciton is formed, e.g., 
from the s— shell conduction band to p— shell conduc- 
tion band 2 ^. Consequently, by applying a lateral volt- 
age larger than the values required to break the symme- 
try, and then using a tt pulse excitation to this state, 
a symmetrically forbidden exciton can be created which 
behaves as the metastable state \m), after lowering the 
applied voltage again. In what follows below, we will as- 
sume that this initial state can be created and focus on 
its evolution after applying STIRAP pulses. 

Initially, the QD is prepared in the metastable state 
\m) and the cavity mode in the vacuum state. An x- 
polarized pump field with a Rabi frequency Q p (t) is ap- 
plied between the metastable state \m) and the biexciton 
state \u). The Hamiltonian of the system, in the rotating 
frame at the pump frequency (interaction picture), can 
be written as 

H = h\\m)(m\ + fiAi|j/)(j/| + ft(Ai + A 2 )\g)(g\ 
+h [fl p (t)\u) (m\ + gi \u) (y\a + g 2 \y) (g\a + H.c] , (1) 

where A p , and Ai (A2) are the detunings of the pump 
field, and the cavity mode with the biexciton (exciton) 
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FIG. 1: (Color online) Schematic diagram for time-bin entan- 
gled photon generation. The transitions from the biexciton 
state |u) to the exciton states \y) and \y) to ground state \g) 
are coupled by a y-polarized cavity mode. The QD is pumped 
from the metastable state \m) to biexciton state through the 
superposition of two input pulses. 



transition frequencies, respectively; H.c refers to Hcrmi- 
tian conjugate. For simulating the dynamics of the sys- 
tem, we perform quantum master equation calculations 
in the density matrix representation. The evolution of 
the QD-cavity system is given by 



p = - \ \h, p\ - \ E \ L l L »p - 2L »P L l + P^M > (2) 



where L M are the Lindblad operators, with terms 
y/ji\m)(u\, y/ji\y)(u\, and yfn\g){y\ corresponding to 
the spontaneous decays, and y/2^\u){u\, \Pfd\vn) (rn\, 
and \P)d\y){y\, corresponding to pure dephasing of the 
biexciton and exciton states. The emission of the photons 
from the cavity mode is given by the Lindblad operator 
y/Ha, where 2k is the decay rate of the leaky cavity. We 
safely neglect the spontaneous decay of the metastable 
state \m) during the evolution, as the lifetime of the 
metastable state is, by definition, very large. 

We numerically solve the optical Bloch equations us- 
ing Eq. ([2]). for density matrix elements ptj — (i\p\j). To 
simplify the nottaion, we use the definitions \Y) = \y, 1), 
\G) = |.g,f), and \G') = \g, 2), where the alpha-numeric 
notation corresponds to the energy state of QD and the 
number corresponds to the cavity photon. The complete 
dynamics of the system is expressed by the following 
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equations of motion: 



Pmrn 
Puu 

Pyy 

PG'G' 
PGG 

Pvv 
Pgg 



PuY 



PYG> 



PmY 



PuG> 



PmG' 



PyG 



—i£lp(t)p um + iflp(t)p m u + llPuu, (3) 

iCl*(t)p um - ifl p (t)p mu + igiPuY - igiPYu (4) 
-^Ixpuu, 

-igiPuY + igiPYu - igi^PG'Y + ig 2 V2pYG> 

-{h + 1 2 )PYY, (5) 
ig 2 V2pG'Y - igzV^pYG 1 - ^KpG'G 1 , (6) 
Wi{PGy + PyG) + 2>KpG'G' + llPYY - UPGG, (7) 
~W2 (PGy ~ PyG) + np Y Y + llPuu - llPyy, (8) 
KPGG + 72 Pyy, (9) 

-iA p p 

urn ^ L pHrnm igiPYm + i^pPuu (10) 



~(7l + —)Pu,n, 

-iA 1 p uY - ittpPmY ~ igiPYY + igiPuu 

+ig 2 v2p u G> - (7i H 2 J* "™' 

-iA 2 pYG> - igiPuG' - ig 2 V2pG>G' 

f 3K + ~f 2 + ld. 



-«32 



V^pyy - ( : 



)PYG r 



-i(A 1 - A p )p mY - iQ*(t)p uY + igiPn 

■ fn i K + 72< 

-ig 2 V2pn 



JmG' 



{id 



mY, 



-i(Ai + A 2 )p u G> - itop(t)pmG' - igiPYG 

+ig 2 V2p uY - {id + K + i 2 )p u G', 
-i{Ai + A 2 - A p )p mG > ~ iCl*(t)p u G' 

+ig 2 V2p mY — (« + -^)PrnG>, 

-iA 2 p yG - ig 2 pGG + ig2P vv 

(k + 72 + 7d) 
o PaG, 



(11) 

(12) 
(13) 
(14) 
(15) 
(16) 

(17) 



When the pump field and the cavity mode satisfy the 
two-photon Raman resonance condition A p wAi + A2, 
the evolution of the population in the QD energy lev- 
els follows the cavity-assisted STIRAP. Here, one must 
also remember that the Stark shifts of energy levels also 
play an important role in two-photon Raman resonance 
condition-^. However, for constant cavity couplings, the 
Stark shifts in energy states remains constant and the 
two-photon Raman resonance condition can be satisfied 
easily by changing the detuning of the pump field (A p ) 
only. 

In Fig. [21 we show the numerical simulations after solv- 
ing Eqns. (|5"))- (fTT)) . The pump field is chosen to be a 
coherent superposition of two time-separated Gaussian 
pulses of the same width, but different amplitudes, which 
can be generated by passing a Gaussian pulse through an 
unbalanced two arm interferometer. The Rabi frequency 
of the pump field is given by Q p (t) = Qi(t) + Q 2 (t — T), 
where fli t2 (t) is the Rabi frequency of each pulse and T 
is the time gap between the pulses. We select a typical 
value of Oi (£) such that the population of the state |m) 
is pumped to the state \y, 1) in STIRAP with probability 




FIG. 2: (Color online) Populations in state \m) (black solid 
line), in state \g) (red dashed line) and in state \u) (blue 
dotted line) during the evolution. Two Gaussian pump pulses 
of different amplitude are employed such that pi = p2 = 
1/2. Using the parameters gi = Q2 = g, k = 2.5g, 72 = 
71 = 10- [, g, ld = l(r 2 5 , Ai = -3g, A 2 = 2g, A p = -1.5g, 
fii(t)/fl = 0.74exp[-(t - 2t p ) 2 /t 2 ], Q 2 (t)/g = 3exp[-(f - 
9.5t p ) 2 /t 2 ], and gr p = 2tt. 



pi = 1/2. Due to the nature of the leaky cavity mode, 
the photon is emitted from the final state \y, 1) and the 
system is evolved into the state \y,0) state. The popula- 
tion in state \y, 0) is transferred to the state \g, 1) through 
the cavity mode. After emitting another photon from the 
state \g, 1), the system finally reaches the ground state 
I g, 0). Thus a photon pair in the early time bin is emit- 
ted during the interaction of the first pulse Oi(t). The 
remaining population in state \m) is similarly pumped 
by another pulse Q, 2 (t — T) and a photon pair is gener- 
ated with probability p 2 — 1 — pi in the late time bin. 
The state of the generated photon pair is thus maximally 
time-bin entangled state. For a QD embedded in a micro- 
cavity, the off-resonant exciton has a spontaneous decay 
rate of the order of 0.1 — 1 ^eV (10~ 3 .g for g = 0.1 meV), 
and the cavity decay condition k » g can be achieved 
easily. We stress that all of these parameters correspond 
closely to those in present day experiment o 17 ' 18 . 

During the time interaction with the pump pulses, the 
population in the upper state \u) remains always less that 
0.1 and the population in \g,2) remains negligible. The 
state of emitted photon pair from the cavity mode can 
therefore be written as 

m)) = [<A(t)4(t) + a\(t - T)4(t - T)J |{0}), (18) 

where |{0}) is vacuum field, (ai(t)) = (a y Y{t)), and 
(02(f)) = {<J g G{t)), for <Jij = |i)0'|- In quantum infor- 
mation protocols, such as entanglement swapping, it is 
essential that the photons in the mode a\ and a 2 should 
not have any other correlation except the time-bin en- 
tanglement. However, in the biexciton-exciton cascade, 
the a\ mode photon is always generated after the emis- 
sion of the a 2 mode photon. Thus the a\ and a 2 modes 
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remain time correlated. This undesirable temporal cor- 
relation becomes negligible for Ti/T 2 >> 1— , where 
is the emission rate of the photon in ai mode. In our 
scheme above, the first photon in mode a\ is generated 
in resonant Raman process, which is emitted with the 
cavity mode decay rate n, and the second photon is gen- 
erated through cavity enhanced spontaneous emission. 
The condition Ti/T 2 >> 1 can therefore be easily sat- 
isfied by choosing a suitably large value of the detuning 
A 2 such that g\nl{n 2 + A|) << k. 



III. TRIPLE COINCIDENT DETECTION OF 
THE PHOTON ENTANGLEMENT, AND THE 
INFLUENCE OF PURE DEPHASING 

Next, we discuss how to measure the entanglement of 
the generated state of the photons. The Concurrence of 
the state ([!%)) is directly related to the coherence of the 
stated. For measuring the degree of entanglement, pho- 
tons from each mode are passed through an unbalanced 
two path interferometer—; the time difference between 
two arms is T, with phase difference <fi, and T is simi- 
lar to the time difference between the two pulses in the 
pump fields. After passing through the interferometers, 
the field operators at the output of the interferometers 
can be expressed as 



The post-selection, for detecting both photon simultane- 
ously after passing through the interferometers, projects 
the state ([TBI) into the state 



a 3 (t) =ai(*) + e*V(t-T), 
a 4 (*) = a 2 {t) + e l *a 2 (t-T). 



(19) 
(20) 



<4(i)4(i) + (1 + e 2l<i >)a\{t - T)a\{t - T) 

+a\(t - 2T)a\{t - 2T)1 |{0}).(21) 



Clearly, the state (|2"Tj) has three terms which are distin- 
guishable in time. The middle term, appearing at t = T, 
provides the information about the entanglement of state 
(fl~8| . For separating different terms in state (|2"Tj). the 
time of detection of photons is measured with reference 
to the pump photons using a triple coincidence detec- 
tion. The probability of triple coincidence detection of 
one photon at the output of each interferometer, and one 
from the input pulse fix, is given by 

/•oc rT bin 

G (3) (t) = / dt' aY'|r!i(0| 2 x 

JO J-Tun 

(al(t' + T)a\{t' + t + T')a 4 (t' + t + T> 3 (i' + r)),(22) 

where Tu n is the width of the time-bins; which is cho- 
sen larger than the biexciton-exciton cascade decay and 
smaller than T. We can simplify the above expres- 
sion for G^\t) using the property of field operators, 
a\{t)a 2 {t — T)|{0}) = 0, as both photons are generated 
almost together in the cascade decay. Subsequently, we 
can simplify the correlation function in (|22j) as 



(a\{t' + r)4(*' + T + T r )a 4 (t' + r + r')a 3 (t' + r)) = (a\(t' + T )a\{t' + r + r')a 2 (i' + r + r>i(f' + r)) 



+ (a[{t' -T+ r)4(t' -T + t + r> 2 (i' -T + t + r')ai(i' - T + r)) 
+ 2cos2(/)(a\(t' + r)4(*' + t + r')a 2 {t' -T + r + r')a x {t' - T + t)), (23) 
I 



which is evaluated for state (fT8|) by applying quantum 
regression formula^. We relegate the details of the 
calculation to the Appendix. 

In Fig. [31 we plot G^(t) for the same parameters 
used in Fig. [21 where the time-dependent populations 
were shown. The computed value of G^ 3 '(r) has three 
peaks centered at r = 0, T, and 2T. The first peak 
at t = correspond to the photons generated in the 
early time-bin that have passed through the short arms 
of the interferometers. Similarly the peak centered at 
t = 2T corresponds to the photon pair generated in the 
late time-bin that have passed through the long arms of 
the interferometers. The central peak at r = T corre- 
sponds to the overlap of the photons generated in early 
time bin and passed through the longer arms in the inter- 
ferometers and the photons generated in the late time- 
bin and passed through the short arms. Thus only the 



central peak contains the information about the entan- 
glement and can be easily selected by choosing a narrow 
time window around r = T. We have also found that the 
required value of T is slightly less than the actual time 
between the pump pulses, which shows that in STIRAP. 
the photons are actually generated before the pump pulse 
reaches its maximum. For the parameters used in Fig. [2l 
the time between pump pulses is lBir/g, but the central 
peak in Fig. [3] is a maximum for T = \Att / g. 

The coherence in the generated state (| can be 
measured by varying the phase <p between the overlap- 
ping amplitudes corresponding to the early and the late 
time bins along the central peak. In FiglH we plot 
the interference pattern produced in the measurement 
of P c = f^T in G (3) (r)dr. The visibility of the interfer- 



ence pattern, defined as V- 
of P c ) / (maximum of P c - 



- (maximum of P c - minimum 

- minimum of P c ), gives the 
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FIG. 3: The triple coincidence correlation of detecting one 
photon at output of each interferometer and one from the 
input pulse fii for T = 14m/g. The other parameters are the 
same as in Fig. 2 
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FIG. 4: The integrated value of the triple coincidence correla- 
tion G' 3 ' (t) along the central peak at r = T. The interference 
pattern appears on changing the phase <j) produced by the in- 
terferometers. 
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FIG. 5: The dependence of the visibility, i.e. entanglement 
of the generated time-bin entangled state of the photon on 
dephasing rate 7d. 



concurrence and the purity of the generated state (fT5|) . 

In Fig. [5j we show the dependence of the visibility 
on the dephasing rate. In the presence of pure dephas- 
ing, the visibility, i.e. concurrence of the generated state 
(Eq. [1~8|) . is strongly inhibited. Further the dephasing in 
the state (Eq. [T8|) only occurs during the interaction of 
the pump pulses. The dephasing during the time gap 
between the pump pulses plays no role as the coherence 
in the state is produced by the coherence between the in- 
put pulses. The inhibition of interference due to dephas- 
ing can be understood as due to pure dephasing (e.g., 
phonon) interactions the information about the photon 
generating in different time bin is partially imprinted in 
the phonon baths and thus as a result of quantum com- 
plementarity the interference is inhibited. However, for 
small dephasing rate jd Rj O.Olgi, the value of visibility is 
larger than 1 / v2, which is required for violation of Bell's 
inequalities^. 

IV. CONCLUSIONS 

We have presented a cavity-QED STIRAP scheme for 
generating a scalable source of time-bin entangled photon 
pairs, and we also investigated the role of pure dephas- 
ing on entanglement. The generated state of the photons 
can be detected by measuring the correlations between 
the pump and the generated photons. We found that for 
small values of pure dephasings, it is possible to achieve 
larger values of entanglement using current working tech- 
nologies. 
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Appendix: Calculation of multi-time correlations 

Here we briefly discuss the method for calculating 
two-time correlation and four-time correlation used in 
Sec. III. We follow the approach discussed by Gardiner 
and Zoller— for evaluating multi-time correlations. The 
required two time correlation can be expressed as 

(a\(t)al(t + T)a 2 (t + r)oi(t)) 

= Tr ia 2 {t + T)ai(t)p(0)a\(t)4(t + r)} (A.l) 

= Tr{a 2 {t + T)p\t)al(t + T)Y (A.2) 

= Tr{a 2 p'(t + T)al} ( A -3) 

where Tr stands for trace and operators a, appear- 
ing without time parenthesis are in Schrodingcr picture, 
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p'(t) = ai(t)p(0)a\(t) = aip(t)a\ is calculated after Finally, the value of four-times correlations used in 
evolving the initial state p(0) = p mm |m)(m| for time t Sec. Ill arc found using Eq. (|A.8[I . 
using Eqs.([3"))- (|17[) and then operating by a x and a\ from 
left and right, respectively. Clearly, p' also follows the 
same equations of motions (f3"T)- (fT7| . Now, using initial 
value p'(t) = aip(t)a\ at time t, and evolving for time r, 
p'(t + t) is calculated. The value of the required correla- 
tion is calculated using Eq. (|A.3|) . A similar straightfor- 
ward approach, considering the times appearing in the a 
operators in ascending order, is applied in evaluating the 
four time correlations as follows: 



where p(0) = p mm \m)(m\, p 1 (t - T) = a x {t - T)p(O) = 
a lP (t - T), p 2 {t -T + t) = a 2 (t -T + r)pi(t — T) = 
a 2 pi(t-T + T), and p 3 (t) = pi(t-T + r)a\(t) = p 2 {t)a\. 
Thus the density matrices p, p±, p 2 , and p% are evolved 
for times, tot-T, t-Ttot-T+r, t-T + T to t, and 
t to t + t respectively. The evolution of p(0) is given by 
Eqs. (|3))- (fl"7|) . while the evolutions of density matrices pi 
for i=l,2,3, follow the similar equations, written for p, as 

f>Yy = "igiPuy ~ ig2\^PG'y + WZpYG 

~{ l2 + K /2 + ld )p Yyi (A.9) 

PG'G = -ig 2 V2pYG + iP2PG'y ~ ^PG'G, (A. 10) 

pGg = -igiPyg ~ ^PGg, (ATI) 
Puy = -iAip U y - iVL p (t)p my - igipYy + ig2PuG 



(a\{t)al{t + r)a 2 {t -T + T )cn{t - T)), 




(A.4) 




PG'y 



-{ 11 + l2 /2 + 3 ld /2) Puy , (A.12) 
iA 2 p G > y - ig 2 V2pYy + igiPG'G 



Pyg 



-{K + y 2 /2 + ld /2)p G , y , (A.13) 

~iA 2 p Y G - igiPuG - ig2"/2~PG'G + ig2PYy 



Pray 



Pyg 



-(K + j 2 /2 + ld /2)p YG , (A.14) 
-iA 2 p yg - ig 2 p Gg - -(72 + ld)p yg , (A.15) 
-i(Ai - A p )p my ~ iQ*(t)p uy + ig 2 p m G 



PmG 



PG'g 



PuG 



PYg 



-(j2/2 + Jd)p my , (A.16) 
— «(Ai + A 2 )p uG - itt p p mG - igipYG + ig2Puy 
~{n/2 + ll + ld ) PuGl (A.17) 
-i(Aj + A 2 - A p )p mG - in*(t)p uG + ig 2 p my 

-\{^ + ld)p mG , (A.18) 
-ig 2 V2p Yg - Kp G , g , (A. 19) 

-iA 2 p Yg - ig 2 V2p G ' g - igip ug 




Pug 



— i(Ax + A 2 )p ug - iQppmg - igiPYg 
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